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DEVELOPMENT AND VERIFICATION OF REAL-TIME, HYBRID COMPUTER 
SIMULATION OF F100-PW-10CX3) TURBOFAN ENGINE 
by John R. Szuch, Kurt Seldner, and David S. Cwynar 
Lewis Research Center 

SUMMARY 

In recent years, there has been increased interest in developing digital, electronic 
controls for airbreathing propulsion systems. Real-time computer simulations of en- 
gines can facilitate the development of these digital controls. The engine simulation 
provides a "test-bed" for evaluating new control laws and for checking and "debugging" 
the control software prior to engine testing. This report describes a real-time, hybrid 
computer simulation of the Pratt & Whitney F100-PW- 100(3) augmented turbofan engine. 
The simulation is intended to support controls research programs involving that engine. 
The simulation has both steady-state and transient calculation capabilities. This report 
describes the modifications that were made to a previously developed simulation of the 
F100-PW- 100(1) engine in order to match the predicted performance of the more ad- 
vanced F100-PW- 100(3) engine. Baseline performance data were obtained from 
Pratt & Whitney's digital simulation of the engine. Data are presented to show that the 
real-time simulation does match the baseline steady-state and transient performance 
over a wide range of flight conditions and power settings. This report also includes the 
equations which describe the F100-PW- 100(3) engine model, FORTRAN listings of the 
digital portion of the simulation, and analog patching diagrams. 


INTRODUCTION 

Over the past several years, aircraft operational requirements have dictated the 
development of gas turbine engines which deliver increased performance over a wider 
operating range. These development efforts have resulted in today's complex, aug- 
mented turbofan engines and will, undoubtedly, lead to increasingly complex, variable- 
cycle engines in the future. 



It is not surprising that, as engines have become more sophisticated, the task of 
controlling those engines to provide safe and stable operation with increased perform- 
ance has also become more difficult. As a result, there has been increased interest in 
applying multivariable (optimal) control theory to the engine control problem (refs. 1 to 
4). These advanced control concepts, however, require the use of a digital computer 
with its inherent precision, logic, and memory capabilities. The digital computer pro- 
vides the control system with more flexibility and versatility then is currently provided 
by hydromechanical controls (ref. 5). 

It has been shown in references 6 to 8 that the use of real-time computer simula- 
tions of engines can facilitate the development of digital controls. The engine simulation 
provides a "test-bed" for evaluating new control laws and for checking and debugging 
of the actual control software prior to engine testing. A real-time simulation also 
allows the control developer to evaluate the timing and sequencing within the digital 
control and to predict the effects of extended digital sampling intervals (ref. 9) on en- 
gine performance. This report describes a real-time, hybrid computer simulation of 
the Pratt & Whitney F100-PW- 100(3) augmented turbofan engine. The simulation has 
both steady- state and transient calculation capabilities and is intended for supporting 
controls research programs involving that engine. The report describes the modifica- 
tions that were made to a previously developed simulation of the F100-PW- 100(1) en- 
gine (ref. 7) in order to match predicted F100-PW- 100(3) engine performance. Com- 
parisons of hybrid simulation and baseline digital simulation data were made over a 
wide range of flight conditions and power settings. The baseline performance data were 
obtained from Pratt & Whitney's digital simulation of the engine. This report includes 
the results of the comparisons, simulation equations, FORTRAN listings, and analog 
patching diagrams . 


ENGINE DESCRIPTION 

The Pratt & Whitney F100-PW- 100(3) engine (fig. 1) is an axial, mixed-flow, aug- 
mented, twin-spool, low-bypass-ratio turbofan. It features improved fan performance 
over the earlier F100-PW- 100(1) version. A single inlet is used for both the fan airflow 
and the engine core airflow. Airflow leaving the fan is separated into two flow streams: 
one stream passing through the engine core and the other stream passing through the 
annular fan duct. The three-stage fan is connected by a through-shaft to the two-stage, 
low-pressure turbine. A ten- stage compressor is connected by a hollow shaft to the 
two-stage, high-pressure turbine. The fan has variable, trailing edge, inlet guide 
vanes. The compressor has a variable inlet guide vane followed by two variable stator 
vanes. Engine airflow bleed is extracted at the compressor exit and discharged through 
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the fan duct during starting. Compressor discharge bleed air is also used to cool the 
high- and low-pressure turbine blades and to power the augmentor turbopump. 

The main combustor consists of an annular diffuser and a chamber with 16 fuel 
nozzles. The engine core and fan duct streams combine in an augmentor and are dis- 
charged through a variable convergent-divergent nozzle. The augmentor consists of a 
diffuser section and five concentric fuel manifolds (zones). 

The engine's bill- of- material (BOM) control system consists of a hydromechanical 
fuel control system and an electronic supervisory control system. The hydromechanical 
fuel control system ( 1 ) meters fuel to the main combustor as a function of the power 
lever angle PLA, the compressor speed N H , the fan discharge total temperature T 13 , 
and the compressor discharge static pressure P g 3 , (2) positions the compressor vanes 
to improve starting and high Mach number characteristics, (3) meters fuel to the five 
augmentor zones as a function of PLA, T 13 , and P g 3 , and (4) controls the nozzle 
area so as to maintain the desired engine airflow during augmented operation. (All 
symbols are defined in appendix A. Numerical subscripts refer to locations in the 
engine (e.g. , fig. 1). ) The electronic supervisory control (1) positions the inlet guide 
vanes to improve inlet distortion tolerance and fan efficiency, ( 2 ) trims the main com- 
bustor fuel flow to satisfy engine limits, and (3) trims the nozzle area to satisfy engine 
airflow requirements . 


ENGINE SIMULATION 
Engine Model 

The mathematical model which described the performance of the F100-PW- 100(1) 
engine was patterned after Pratt & Whitney's digital simulation (CCD 1015) of that en- 
gine and was reported in reference 7. Subsequent modifications were made to elements 
of that model to match the performance of the F100-PW- 100(3) engine as predicted by 
the corresponding digital simulation (CCD 1103-1.0). Those modifications are de- 
scribed in the following section. The basic structure of the mathematical model was 
not changed, however. Figure 2 contains a computational flow diagram of the F100-PW- 
100(3) real-time simulation. Appendix B contains a complete list of equations which 
define the simulation model. Table I contains a list of engine design parameters for 
the FI 00- PW- 100(3) simulation. 


Simulation Modifications 

The equations describing the mathematical model of the F100-PW- 100(3) engine 
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were implemented on the Lewis Research Center’s hybrid computing system. This 
system consists of an EAI model 640 digital computer, a model 680 analog computer, 
and a model 681 analog computer. The split of the computational load between the 
digital and analog computers was basically the same as that employed in the earlier 
F100-PW- 100(1) simulation (ref. 7). The modifications that were made to the digital 
portion of that simulation included the following. 

First, in the earlier F100-PW- 100(1) simulation, all analog inputs to the digital 
computer were sampled at the beginning of the digital cycle, and all outputs to the 
analog computer were transferred after all of the digital calculations were completed. 
From a dynamic standpoint, this proved to be the worst approach since it resulted in the 
greatest effective time delay - hence, phase shift (refs. 10 to 12). In the F100-PW- 
100(3) simulation, the analog inputs to the digital computer are sampled as needed, and 
the resultant digital data are transferred to the analog computer as soon as they are 
available (on a component by component basis). This approach results in a significant 
reduction in the phase shift associated with individual computational loops (ref. 12) 
since the calculation time for each loop is much less than the total update time. Auxil- 
iary calculations such as the calculation of engine thrust and surge margins contribute 
only the total update time. 

Second, the fan and compressor performance maps represented by equations (Bl), 
(B2), and (B15) were based on axial vane positions in the earlier F100-PW- 100(1) simu- 
lation. This necessitated shifting the map data when operating at low corrected speeds 
where the vanes are cambered. To minimize this shifting in the F100-PW- 100(3) 
simulation, the fan and compressor maps were regenerated with the vanes on their 
nominal schedules. Therefore, no shifting of the map data is required when the vanes 
are on their schedules. During transients and other off-schedule conditions the required 
shifting of corrected airflows is accomplished by equations (B3) and (B16) with the 
shifts computed from bivariate functions of corrected speed and vane position. 

Third, to better match the F100-PW- 100(3) baseline digital data over the entire 
flight envelope, an empirical Reynolds number effect on fan performance was added. A 
shift in the fan corrected airflow (eq. (B3)) is computed as a piecewise linear function 
of the Reynold's number index (eqs. (B4) and (B5)). 

Fourth, all the bivariate component performance maps and shift functions were re- 
generated to match the predicted FI 00- PW- 100(3) steady-state performance. Those 
curves are shown in figures 3 to 8. 

Fifth, surge margin calculations were added for the fan (eqs. (B9) to (B12)) and the 
compressor (eqs. (B19) to (B22)). In each case, the critical pressure ratio was fit by 
a quadratic function of the corrected airflow at each of the extreme vane positions. The 
quadratic functions were based on fits of digital simulation data. 

Sixth, the exhaust nozzle exit area A g was fit by a linear function (eq. (B52)) of 
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the nozzle throat area A ^ and the nozzle inlet temperature for each of two ranges 
of flight Mach number. In the earlier F100-PW- 100(1) simulation, the effect of nozzle 
heating (T^) was not considered. 

Seventh, the inlet calculations of P 2 and T 2 were eliminated in the FI 00- PW- 
100(3) simulation. These variables are transferred as input to the digital portion of the 
hybrid computer from the analog computer, thus allowing operation of the simulation 
during changes in the flight condition. 

Lastly, in the earlier F100-PW- 100(1) simulation, the same digital program that 
was used to perform the required calculations was also used for input and scaling of 
component performance data and for setup of the analog consoles. In the F100-PW- 
100(3) simulation, these functions are performed by separate digital programs. The 
scaled, component performance data are shared by the data input program, the main 
digital program, and the function generation routines through the use of COMMON 
blocks . 

Appendix C contains a FORTRAN listing of the digital portion of the F100-PW- 100(3) 
real-time hybrid computer simulation. Reference 6 contains a detailed discussion of 
the digital program structure including the MAP2 and MAP2L function generation rou- 
tines. 

Modifications were also made to the analog portion of the F100-PW- 100(1) real-time 
simulation. These included the following: first, to better match the predicted F100-PW- 
100(3) augmentor pressure drop, the pressure drop was computed using the total aug- 
mentor gas flow (including augmentor fuel flow) and the discharge temperature T^ 

(eq. (B41)). In the earlier F100-PW- 100(1) simulation, the effects of augmentor fuel 
flow and its associated energy release were not included. Second, the augmentor effi- 
ciency and duct pressure drop curves (eqs. (B43) and (B45)) were regenerated to better 
match digital simulation data over the entire flight envelope. The new curves are 
shown in figures 9 and 10. Third, the exhaust nozzle discharge coefficient was fit by a 
piecewise- linear function of the nozzle pressure ratio (eq. (B49)) having more segments 
than in the earlier FI 00- PW- 100(1) simulation. Fourth, in the F100-PW- 100(1) simula- 
tion, the fen discharge (core side) and compressor discharge specific heats were as- 
sumed to be linear functions of the corresponding temperatures for the purpose of com- 
puting the required torques. The intercepts of the linear functions were adjusted at the 
military power setting to match the baseline rotor speeds at each flight condition. 

Changes in flight condition could not realistically be accomplished, however, because of 
the wide range of intercept values. To eliminate this problem in the FI 00- PW- 100(3) 
simulation, the intercepts were also fit by linear functions of the fan inlet pressure and 
temperature (eqs. (B80) and (B85)). Only slight adjustments of the resultant intercepts 
were then required to match rotor speeds after changes in the flight condition. Realis- 
tic changes in the flight condition could be accomplished with fixed intercept values. 
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Lastly, calculations of the fan discharge (P - P )/P, for both the duct and core sides, 
were added to the real-time simulation (eqs. (B87) and (B88)). Figure 11 shows the 
functional relation between (P - P_)/P and the compressible flow parameter. 

Appendix D contains the analog patching diagrams for the F100-PW- 100(3) real- 
time, hybrid computer simulation. 


Simulation Requirements 

The digital portion of the F100-PW- 100(3) real-time simulation consumed 12 440 
words of core storage (including data). The supplemental data input program consumed 
7144 words of core storage (including data). The digital computer update time, which 
was approximately 7. 5 milliseconds, resulted in stable, real-time operation. 

Both analog computers were fully utilized. For example, the full complement of 
24 multipliers on the 680 analog computer and 30 multipliers on the 681 analog computer 
were used. In addition, the full complement of 6 digital to analog multipliers on the 
680 analog computer was used. A total of 189 potentiometers was required. The eight 
digitally set, univariate function generators available on the 681 analog computer were 
also used. 


RESULTS AND DISCUSSION 

The usefulness of the F 100- PW- 100(3) real-time, hybrid computer simulation de- 
pends on its ability to accurately represent the physical engine over the desired range 
of operation. It should match the steady-state and the transient engine performance 
for power settings from idle to full augmentation (maximum thrust) at altitudes and 
flight speeds within the engine operating envelope. Figure 12 shows the flight conditions 
selected for evaluation. 


Procedure 

As previously stated, the basis for comparison of hybrid simulation data was the 
engine manufacturer's digital simulation of the engine (CCD 1103-1. 0). That simula- 
tion also included a simulation of the BOM control logic. The digital simulation of the 
engine and control could be run in either a steady-state (fixed PLA) or transient (time- 
varying PLA) mode. 

All of the selected flight conditions (fig. 12) were first run in the steady-state mode 
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with PLA ranging from the minimum allowable setting to the maximum thrust setting of 
130°. The minimum allowable setting (idle) was dictated by the control and was based 
on inlet airflow requirements and minimum combustor pressure limits. For the three 
subsonic conditions having altitudes lower than 10 kilometers, the idle setting was 20°. 
For the 13. 72 km/k^ = 0. 9 condition, the idle setting was 30°. For the three super- 
sonic conditions, power settings lower than 83° were not permitted. 

The steady-state data obtained from the digital simulation included (1) values for 
the control variables such as main combustor fuel flow and (2) values for selected en- 
gine variables such as fan speed. The hybrid computer simulation was then evaluated 
in steady state at each flight condition and power setting by setting the analog control 
inputs at the appropriate values and then recording the resulting values of the selected 
engine variables. In this way, the simulation could be evaluated without requiring a 
separate control simulation. This open-loop approach ensured that observed differences 
between hybrid and digital simulation data were attributable to the hybrid simulation and 
not to control simulation errors. The following section compares the F100-PW- 100(3) 
hybrid simulation and baseline digital steady-state data at the selected flight conditions. 

The four subsonic flight conditions (fig. 12) were also run on the digital simulation 
in the transient mode. In each case, the PLA was initialized as its minimum value and 
then stepped to 83° (at t = 0 sec). The 83° setting was maintained for 10 seconds at 
which time the PLA was stepped down to its minimum value. The transient data obtained 
from the digital simulation provided time histories of both the control variables and the 
selected engine variables. The hybrid simulation was then evaluated for transient 
operation by scheduling the control inputs to the hybrid simulation to match the digital 
time histories. The resulting engine response data were recorded and subsequently 
compared with the digital results. As in the steady-state evaluation, this open-loop 
approach eliminated the need for a separate, real-time control simulation, and it 
allowed the isolation of engine simulation errors. A following section (p. 9) compares 
the transient data obtained with the F100-PW- 100(3) hybrid computer simulation with the 
corresponding baseline digital data. 


Steady- State Simulation Results 

The verification of the steady-state performance of the FI 00- PW- 100(3) real-time 
simulation was accomplished by operating the simulation in an open- loop manner at each 
of the flight conditions shown in figure 12. At each power setting, the values of the 
main combustor fuel flow, exhaust nozzle area, fan inlet guide vane angle, compressor 
stator vane angle, and augmentor fuel flow were set to match the baseline digital values. 
The engine variables selected for the steady- state comparison were fan speed, com- 
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pressor speed, main combustor pressure and temperature, net thrust, fan-tip pressure 
ratio, total fan corrected airflow, compressor pressure ratio, and compressor cor- 
rected airflow. Agreement of hybrid and baseline digital values for these variables 
would represent good, overall steady-state verification of the real-time simulation. 

Figures 13 to 17 contain plots of the hybrid and digital steady-state data at the 
selected flight conditions. For convenience, the engine variables were plotted against 
the PLA which corresponded to the set of control variables. It should be noted that, at 
each flight condition, the plot scales were expanded to match the observed range of the 
data. The scale expansion was most significant at the supersonic flight conditions. 

Prior to recording steady-state, hybrid simulation data at each flight condition, the 
fan and compressor discharge specific heats were adjusted to achieve a match of base- 
line rotor speeds at the 83° power setting as shown in figures 13 and 14. Agreement of 
hybrid and digital simulation data at other power settings, however, was dependent on 
the accuracy of the individual component models (i.e. , fan, compressor, nozzle, etc. ). 
Agreement of hybrid and digital data at these conditions would serve to substantiate the 
simulation simplifications that were required to achieve real-time operation. 

Figure 13(a) shows excellent agreement of fan speed along the entire sea-level/ 
static operating line. This was attributed to the fact that the temperature- sensitive 
specific heat relations (eqs. (B80) and (B85 )) were established at the sea- level/static 
condition. Figures 13(b), (c), and (d) show good agreement at the other flight conditions 
although some discrepancies were observed in the midpower range at the higher alti- 
tudes. These errors were less than 2 to 3 percent of the design fan speed. Figures 
13(e), (f), and (g) show the results of the fan speed comparison for the 6. 096 km/M n = 

1. 8, 12. 19 km/M n = 2.2, and 17. 83 km/M n = 2. 15 conditions, respectively. For these 
conditions, the comparison was limited to power settings of 83° and above. For these 
conditions, good agreement in fan speed was also observed. The errors were generally 
less than 2 percent of the design speed and were attributed to the assumption of constant 
gas properties in the hybrid simulation model of the exhaust nozzle. 

Figure 14 compares the hybrid and baseline digital results for the compressor 
speed. For all power settings below 83°, the observed errors were less than 3. 5 per- 
cent of the design speed. For the supersonic, augmented operating points, the errors 
were less than 1. 1 percent. 

Figures 15 and 16 compare the baseline digital and hybrid simulation values for the 
main combustor and temperature, respectively. As in the case of the rotor speeds, ex- 
cellent agreement was observed for power settings below 83° throughout the operating 
envelope. Errors were generally less than 3. 5 percent of the design value. For the 
supersonic, augmented conditions, errors in the main combustor pressure and temper- 
ature were less than 2. 1 percent. 
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Figure 17 shows the comparison of net thrust for the selected flight conditions. In 
the thrust calculation, a constant velocity coefficient was assumed and, as in the calcu- 
lation of the exhaust nozzle flow, constant gas properties were assumed. Even with 
these simplifications, generally good agreement between hybrid and baseline digital 
values for thrust was observed. Excellent agreement was observed at the subsonic con- 
ditions. However, differences of up to 9 percent of the design maximum thrust were 
observed at supersonic, augmented conditions (see fig. 17(f)). 

The fen and compressor operating lines are shown in figures 18 and 19, respec- 
tively. Good agreement between the baseline digital and hybrid simulation data was ob- 
tained for all flight conditions. Some discrepancy in fan corrected airflow (about 
2. 5 percent of the design value) was observed at the 13. 72 km/M n = 0. 9 condition for 
high power settings. This error is attributed to Reynolds number effects, since this 
condition represented a lower Reynolds number index than the other selected subsonic 
conditions. An attempt was made to incorporate a fan airflow shift as a function of the 
Reynolds number index, but its adequacy was limited by a restriction on the maximum 
digital update time allowable for real-time operation. Good agreement was also ob- 
served for the compressor operating line. A maximum error of 4 percent in corrected 
airflow and 3.5 percent in pressure ratio was observed at the 3. 048 km/M n = 0. 9 con- 
dition at the idle power setting. 

The results presented in figures 13 to 19 indicate that the hybrid simulation ade- 
quately matches the baseline digital simulation in representing the steady-state behavior 
of the FI 00- PW- 100(3) engine. The hybrid simulation errors (relative to the digital 
simulation) that were observed were sufficiently small so as to indicate that the hybrid 
simulation could be used to evaluate steady-state control functions such as speed regu- 
lation, temperature limiting, and surge protection. 


Transient Simulation Results 

The previous section demonstrated the capability of the F100-PW-100(3) real-time, 
hybrid computer simulation to predict the steady-state performance of the engine. The 
hybrid simulation must also predict the transient performance of the engine so as to 
serve as a tool for developing research control systems. The comparison between 
baseline digital and hybrid simulation data is presented for the four subsonic flight con- 
ditions shown in figure 12. The subsonic conditions were selected since they permitted 
variations in the PLA below the 83° setting. 

The five control inputs to the hybrid simulation were scheduled as functions of time 
to match baseline digital values for a power lever ramp (slam) from the idle setting to 
the 83° power setting. The schedules included a power lever cutback (chop) from 83° 
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to idle 10 seconds after the initiation of the transient. The open- loop operation was 
selected for the transient evaluation to allow the isolation of simulation errors from 
potential control simulation errors. 

Figure 20 shows the comparison of baseline digital and hybrid simulation responses 
to the simulated power lever movement at the sea-level/static condition. The responses 
of fan speed, compressor speed, compressor discharge pressure, main combustor tem- 
perature, and thrust are presented. A slightly higher fan speed overshoot (about 
0.8 percent) and faster deceleration were observed for the hybrid simulation (fig. 20(a)). 
The hybrid simulation response of compressor speed (fig. 20(b)) was slightly faster 
for both acceleration and deceleration. This was attributed to the simplifications used 
in modeling the compressor temperature ratio (torque). Figure 20(c) shows the re- 
sponses of the baseline digital and hybrid simulation values of compressor discharge 
pressure. The responses match quite well except for a discrepancy (about 3. 5 percent) 
at the end of the acceleration. This error was attributed to the simplified compressor 
temperature ratio calculation in the hybrid simulation. Figure 20(d) compares the digi- 
tal and hybrid simulation responses of the main combustor temperature. The most 
notable difference in the responses was that the hybrid simulation resulted in a 4 percent 
lower temperature rise at the start of the acceleration. This temperature difference 
was maintained throughout the acceleration. The hybrid simulation response exhibited 
no temperature overshoot while the baseline digital response overshot the final temper- 
ature by 3 percent. These differences could have been caused by any number of sim- 
plifying assumptions in the hybrid simulation. One of these was the absence of any 
pressure effects in the main combustor efficiency calculation. Figure 20(e) shows a 
similar discrepancy (about 2 percent) in the thrust responses at the end of the acceler- 
ation. The observed discontinuity in the hybrid simulation thrust was due to a simpli- 
fication in the thrust calculation. Because of limits on the digital calculation time, it 
was not possible to accurately model the nozzle performance when normal shocks would 
exist in the nozzle's divergent section. Since this condition only exists at low altitude, 
low speed, low power conditions, it was not considered to be a serious problem. The 
assumption was made that flow at the nozzle throat would be either subsonic or sonic 
(with the shock expelled). Therefore, the discontinuity represented a switch from sonic 
to subsonic flow or vice versa. 

Figures 21 to 23 show comparisons of baseline digital and hybrid simulation tran- 
sients for the 3.048 km/k^ = 0. 9, 9. 144 km/M n = 0. 9, and 13. 72 km/M n = 0. 9 condi- 
tions. In general, the transient results at these conditions were similar to the results 
obtained at the sea-level static condition. That is, the hybrid responses exhibited more 
fan speed overshoots (figs. 21(a), 22(a), and 23(a)) and slightly faster compressor speed 
responses (figs. 21(b), 22(b), and 23(b)). Figures 21 to 23 also reflect some of the 
steady-state differences that were discussed in the previous section. Examples of this 
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are lower main combustor temperature and higher net thrust at 83° PLA for the 
13. 72 km/M n = 0. 9 condition (figs. 23(d) and (e)). 


SUMMARY OF RESULTS 

An existing real-time, hybrid computer simulation of the Pratt & Whitney F100- 
PW- 100(1) turbofan engine was modified to match the predicted performance of the 
F100-PW- 100(3) turbofan engine. The basis for the simulation modifications was the 
engine manufacturer's digital simulation (CCD 1103-1.0) of the F100-PW- 100(3) engine. 
The resulting hybrid computer simulation was implemented on the Lewis Research 
Center's EAI model 640 digital computer, model 680 analog computer, and model 681 
analog computer. The digital computer update time was approximately 7. 5 milliseconds 
and resulted in stable, real-time operation. The digital portion of the simulation re- 
quired 12 440 words of core storage (including data). Both analog computers were fully 
utilized. 

The real-time, hybrid computer simulation of the F100-PW- 100(3) turbofan engine 
was evaluated at a number of subsonic and supersonic flight conditions. The evaluation 
covered both steady- state and transient operation. The resulting hybrid simulation data 
were compared with baseline digital simulation results. 

The steady-state evaluation showed that the hybrid simulation generally matched 
the baseline digital simulation within 4 percent over the F100 flight envelope. Better 
agreement was noted at the low altitude/low Mach number conditions since the hybrid 
simulation was designed to match sea- level/static data from the digital simulation. 

The transient evaluation covered large changes in the pilot command at subsonic 
flight conditions. In general, the agreement between hybrid and digital results was 
good. The hybrid simulation did exhibit slightly more fan speed overshoot during ac- 
celerations. Also, the response of the hybrid-simulated compressor speed was faster 
than the digital response. Some of the observed transient differences could be attributed 
to 2 to 4 percent errors in the steady-state values at the initial, idle power settings. 

The results of the evaluation indicated that the real-time, hybrid computer simula- 
tion of the F100-PW- 100(3) turbofan is suitable for use in the development and evalua- 
tion of digital control systems. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 28, 1977, 

505-05. 
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APPENDIX A 


SYMBOLS 


A 

C d 

C P 

F 

FN7 

f i 

f/a 

GVIPOS 

Sc 

HVF 

HVSPOS 

hp 

Ah 

I 

J 

k ab 

k b 

k blwht 

k blwlt 

K i 

% 

K PR5 

l 

M n 

N 

P 

P/P 


cross-sectional area, cm^ 
nozzle flow coefficient 

specific heat at constant pressure, JAg-K 
thrust, N 

nozzle flow function 
functional relation, i = 1 to 14 
local fuel-air ratio 
fan inlet guide vane position, deg 

2 

gravitational conversion factor, 100 cm-kg/N-sec 

heating value of fuel, J/kg 

compressor stator vane position, deg 

turbine map enthalpy drop parameter, JAg-K ' -rpm 

turbine enthalpy drop, JAg 

2 

polar moment of inertia, N-cm-sec 

mechanical equivalent of heat, 100 N-cm/J 

2 2 /■ 2 4 

augmentor pressure loss coefficient, N -sec /kg -cm -K 

main- combustor pressure loss coefficient, N -sec Ag -cm -K 

fraction of high-pressure-turbine cooling bleed that is performing work 

fraction of low-pressure-turbine cooling bleed that is performing work 

component temperature rise coefficient, i = 1 to 16 

1 /2 

nozzle flow constant, kg-K / /N-sec 

low-pressure-turbine discharge pressure loss coefficient 

length, cm 

Mach number 

rotational speed, rpm 

2 

total pressure, N/cm 
pressure ratio 
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PLA power lever angle, deg 

2 

P_ static pressure, N/cm 

s 

Q torque, N-cm 

R^ gas constant of air, 2. 8699x10^ N-cm/kg-K 

REI Reynolds number index 

SMC compressor surge margin 

SMF fan surge margin 

T total temperature, K 

T/T temperature ratio 

t time, sec 

3 

V volume, cm 

W stored mass, kg 

w mass flow rate, kg/sec 

w c corrected mass flow rate, kg/sec 

w turbine map flow parameter, kg-K-cm /N-rpm-sec 

r 

y specific heat ratio 

6 total pressure relative to sea- level conditions 

r) efficiency 

6 total temperature relative to standard- day conditions 

t time constant, sec 

Subscripts: 

AB augmentor 

ax axial vanes 

B main combustor 

BLC customer bleed 

BLHT high-pressure- tur bine cooling bleed 

BLLT low-pressure-turbine cooling bleed 

C compressor 

cm cambered vanes 

cr critical 



D fan duct 

des design 

e nozzle exit plane 

F fuel 

FAN fan 

H high 

HT high-pressure turbine 

I inlet 

ID fan hub (core) 

i initial conditions 

j engine station (fig. 1); j = 0, 2, 2.1, 2.2, 3, 4, 4.1, 5, 6, 7, 8, 13, 16 

j' entrance to volume at station j; j = 3, 4, 4. 1, 6, 7, 13, 16 

L low 

LT low-pressure turbine 

M map 

m measured 

N nozzle 

n net 

OD fan tip (bypass) 

SUB subsonic 

SUP supersonic 

TPBL turbopump bleed 
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HI H 


APPENDIX B 


SUMMARY OF EQUATIONS 


(*c). 


= ft 


' P 13 N L 


FAN, M *yP 2 ^ 


= fo 


' P 13 N L 


'FAN, ID 


2 ie 


Wo = < 


(* C ) 


FAN, M 


N, 


1. + f 12 — GVIPOS 


fl 3 (RED 


REI = 


5 2 (T 2 + 110.33) 
398.5O(0 2 ) 2 


f 13 - °* 0 


if REI > 0. 61168 


'N 


1. 0069 - 1. 6461 REI if 0. 28404 < REI < 0. 61168> 


= 3. 2240 - 9. 4517 REI otherwise 


P 2. 1 - P 2 . 2 


'FAN, ID 


FAN, OD 


13 




K, — ^ + K„ if-^> 2.851 


Pr 


13 


= Ko — — + K if 1.803 < — ^ < 2.851 
P 2 P 2 


^13 

= Kg + Kg otherwise 

P 2 




(Bl) 


(B2) 


(B3) 


(B4) 


(B5) 


(B6) 


(B7) 
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(B8) 




FAN, OD 


SMF = 


13 


xr, FAN 2 J 


cr, FAN 


2 5. + GVIPOS \ /P\ _ /GVIPQS ) /P\ 

/ cr, FAN v 25 ' ' ' P 'cr,ax, FAN ' 25 ‘ ' 'Hr, cm, FAN 


/p\ 

W 


. 9 

„ W 9 0 9 4 W 9 V0 9 

2. 7372x10" 4 —^+ 5. 1591xl0 -4 — — — - + 0.77163 


cr, ax, FAN 


• 2 a 

a W 9 tf 9 n w 9 W d 9 

= 5. 3717x10 4 — — - 3. 0788x10“ ** T + 1. 7258 


Wr, \ 6 


cr, cm, FAN 


-) ■ K 7 (I) + K 

'FAN, ID '^'FAN.ID 


= K r 


^ * 


P/ 


K 


FAN, ID 


10 


~\ 


if ^ 


> 2.889 


'FAN, ID 


if 1.82 < (-1 < 2.889)* 

^FAN, ID 


— j + K 12 otherwise 

' P 'FAN,ID 


J 


T 2. 1 “ T 2 


*■© 


FAN, ID 




— f r 


P 3 N H 


C ’ M ^ P 2-2 


(B9) 


(BIO) 


(BU) 


(B12) 


(B13) 


(B14) 


(B15) 
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w< 


1.2 '(»c), 


C, M 


/ N H \ 

1. + fJ - 5 , HVSPOSl 

W 0 2.2 / 


2.2 




(B16) 


“) = K 13^- +K 14 

A/ C 2.2 2.2 


Kjg + Kjg otherwise 


2.2 


(B17) 


t 3'-^; c ‘2.2 


(B18) 


SMC = 


/ p\ 

V 

> P ir,C 

P 2. 2 


cr, C 


(B19) 


'cr,C 


= /40 + hvspqs\ /p\ _ /: 

V 44 A p 4r,ax,c V 


HV SPQS - 4 \ /P 
44 / \P/ 


cr, cm, C 


(B20) 


• 2 

Wn ftfin n W, 

= 0. 00972 + 0. 65235 


k 2.2^2.\ 


- 11.7048 


cr, ax, C 


2.2 


2.2 


(B21) 


= 0. 11524 


cr, cm, C 


w 2 . 2 0 2.2 


2.2 


2. 9650 W2 - 2 ^ 2 - 2 . + 24. 779 

6 2.2 


(B22) 


w 


w 


BLHT = ° - 01621 w 2. 2 

(B23) 

BLLT = 01436 w 2 . 2 

(B24) 

W TPBL = 0 

(B25) 
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(B26) 


(w p ) - f. 

\ p/jjp 


P 4. 1 N H 


5 1 p ’ 

M t/T 


*4 - (*p) 


P 4 N H 


HT T, 


(hp)jjrp - fg 


P 4. 1 N H 




(^h)jjT - (hp)jjrp 1^4 Njj 


(wT) 


4. V 


w. 


c p, 4 T 4 

(Ah)jjrp 

+ W BLHT 

c p, 3 T 3 K BLWHT^ Ah ^HT 

c p,4.1 

o 

y 

1 

c p,4.1 C p,4.1 


(*p) 


= fr 


P 5 N L 


LT 


’ 4 - 4 fiZl, 


*4. 1 * (*p) 


P 4. 1 N L 


LT T 


4.1 


(hp) LT - fg 


P 5 N L 


' P4 - 1 fiZy 

(Ah) LT = (hp) LT y T 4 4 Nj 


(wT) fi , = Wl3Cp>16Tl6 + w 


"P-6 


4. 1 


c p, 4. 1 T 4. 1 _ ^ Ah ^LT 


'P> 6 


“P,6 J 


+ w 


BLLT 


c p, 3 T 3 _ K BLWLT^ Ah ^LT 


L c p,6 


~P> 6 


(B27) 

(B28) 

(B29) 

(B30) 

(B31) 

(B32) 

(B33) 

(B34) 

(B35) 
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(B36) 



P 5 “ K PR5 P 6 



(wT) 7 , _-Pi^ w 6 T 6 + ^7 a B W F, 7 

C P,7 C P,7 

^AB = f 9 





c p, 16 T 16 - c p, 13 T 13 



(B37) 

(B38) 

(B39) 

(B40) 

(B41) 

(B42) 

(B43) 

(B44) 

(B45) 

(B46) 

(B47) 
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FN7 = 0. 2588 


if/P) <0.53 

K 


C"V-©: 


0.2857 


otherwise 


C d 7 = 0.97031 


if(^ < 0.3334 

H 


1. 0645 - 0.28248(- 


N 


-(II 


if 0.3334 <(-l < 0.3851 

\ 13 , 

'N 


0.98563 - 0.077661- 
VP 


N 


if 0.3851 <[-) < 0.4783 V 

H 


0. 89488 + 0. 11209 if 0. 4783 <0.8184 


'N 


0. 77344 + 0.26048/ — ] otherwise 

V p / N 


N 


w? =K N P 7 (FN7)C d}7 


[a, - t n (T 7 )] 


* 11 = 0 


if T ? < 800 


= 0. 46870 T ? - 374. 97 if 800 == T ? < 860. 6 

= 0. 12681 T ? - 80. 710 if 860. 6 < T ? < 1200 

= -0. 03577 T ? + 114. 39 if 1200 < T 7 < 1993 


= -0. 13355 T 7 + 309. 23 otherwise j 


(B48) 


(B49) 


(B50) 


(B51) 
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Ag = 1.4693 A ? + 0. 19333 T ? - 1048.2 if M n < 1. f| 

= 1.6175 A rj + 0. 12008 T 7 - 753. 16 otherwise 

v2 


(I) 


A /A 

= 4. 7317— - 1.64861 — | - 2.5089 
SUB A 7 \ A 7; 


(I), 


An /An ' 

= 4.7317 — + 1.64861— 1 + 3.5655 


SUP 
P 


s^e _/P 


P 7 V P 'N 


a “(lug 


© 

V ^SUP 


N ' 'SUB 


otherwise 


F g = 46. 405 w ? I T ? ( 0. 2578 - 0. 2578 

P 7 


if l*\ 


P/ N ' P; SUB| 



F „ = I V 20 - 041 ’ i '2 M nl / To 


W, 


-Jg *"2.2 " "BLHT - "BLLT - "TPBL - *BLC ‘ *3> * + W 3. i 


T 3=7X‘* T 3’- T 3> dt + T 3,i 


(B52) 

(B53) 

(B54) 

(B55) 

(B56) 

(B57) 

(B58) 

(B59) 
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(B62) 


P 4=^/ 0 V T) 4'- i 4 T 4] dt+P 4,l 


4 


T _ V 4 P 4 
4 R A W 4 

(B63) 

W 4. 1 = f Q (v} 4 + *BLHT " *4. 1> + W 4. 1, i 

(B64) 

1 = t 4A fn [ (wT) 4.r‘ W 4.1 T 4.l] dt+P 4.1,i 
v 4. 1 u 

(B65) 

rp V 4. 1 P 4. 1 

R A W 4.1 

(B66) 

W 13 =^ t ^2-^2.2-^13 )dt + W 13,i 

(B67) 

T 13 = T~f n ^ T 13 t " T 13 ^ dt + T 13,i 
t 13 U 

(B68) 

p 13 = V A W 13 T I3 
v 13 

(B69) 

C (^13 + *4. 1 + *BLLT + ™TPBL " V dt + W 6, i 

(B70) 



(B71) 


p 6 =^r [<* T V - W 6 T e] dt + P 6,i 

v 6 u 


T 6 


V 6 P 6 

R A W 6 


Wr 


=L (w 6 + w F>7 - w 7 )dt + W 7>i 


P 7 


R A y 7 


^ [(wT) ?t - w 7 t 7 ] dt + p 7>i 


T 7 


V 7 P 7 

r a w 7 


w 


13 



'Ag c \ r t 


l / 'O 
D 


(P 16' " P 16^ dt + w 13,i 


W c 



f\ P ? , - P 7 )dt + w 6>i 


AB 


(NQ)^ (Ah) HT (w 4 + k blwHT^BLHT^ 


(NQ) C =— c p,2.2 w 2.2 

71 


"P> 3* 

|_\ C p,2.2, 


l T 3’ ' T 2. 2 


"P> 3* \ 


"P, 2.2/ 


= 1. 0815 + 0. 00008 T 3 - 0. 00033 T 2 - 0. 00038 P 2 


* Adjusted, if necessary, to match rotor speeds at PLA = 83°. 


(B72) 

(B73) 

(B74) 

(B75) 

(B76) 

(B77) 

(B78) 

(B79) 

(B80) 
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(B81) 


N 


H 


3o r 
"Wo 


(NQ)jjrp - (NQ)^-. 


N 


H 


dt + N H,i 


( n q) lt ( Ah )LT^ w 4. 1 + K BLWLT w BLLT^ 


30J 


C n 9^9 ~ ™9 o) [ P ? ^ T 1S , - Tj 


(NQ)fAN, OD = — u p, 2^ w 2 - w 2. 2 > l"T 1 13’ " A 2) 

P,2 


(NQ) F AN, id " ~ C p, 2. 2 W 2.2 


P> 2 . 2 I rp rp 

' a 2.2 1 2 


L\ C P,2 


j C P;2.2 \ _ 1 0515* _ 0.00011 T 2 2 + 0.00012 T 2 - 0.00159 P 2 
\ C P,2/ 


n l = — — r j 

'Wo L 


( n q) lt - ( n Q)fan, op ~ ^ nq ^fan, id 


n t 


dt + N t < 

Li, X 


^ P 13 ~ P s, 13\ _ f / w 13 /^13 


14 


13 


P 13 A 13 


2.1 ~ ^s, 2.1 ^ _ f f w 2. 2 V^2. 2 


fPo T - P. 

P 2.1 


l 14| 


P 2. 2 A 2. 2 


* Adjusted, if necessary, to match rotor speeds at PLA = 83°. 



(B82) 

(B83) 

(B84) 

(B85) 

(B86) 

(B87) 

(B88) 
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APPENDIX C 


DIGITAL PROGRAM 
FORTRAN Listing 


C****+HDC VAR I ABLES 

SCALED FRACTION X0, XI, X2, X3, X4, X5, X6, X7, XS, X3, X10, Xll 
SCALED FRACTION XI 2, X13, X14, X15, X16, X17, X18, X19, X20, X21, X22 
SCALED FRACTION X23 
C*****DAC VAR I ABL.ES 

SCALED FRACTION VO, Vi, V2, V3, V4, V5, VS, V7, VS, Y9, ViO, Vll, V12 

Ct +-+++ other variables 

SCALED FRACT I ON XVALS < 10, 8,6), WALS < S, 6 > , ZVALS < 10, 8, 12 ) , V I < 24 ) , 

1 V9, VIO, SSQRT, MAF'2, MAP2L, V3, V4, V5, AR, F'RSUB, DF'R, PE, PRE, 

2 F ' RSUP ' FF:D ' VS, V7, VS, RT4, RT41, AE, A, DV7, Vll, REV, PRF1, PRF2, PRC1, 

3 PRC2, VSSQ, V7SQ 

COMMQN/MAF'S/XVALS, WALS, ZVALS, I X < 6 ) , JV<6), NX '16), NY<6>, KX < 6 ) , KV< 
COMMON/ I DAC/V I 


LOGICAL SENSH, READY 
REAL MO 

CALL QSHY I N v I ERR, 680, 630 > 

TVF'E 5 

5 FORMAT <3X, 21HTVPE DATE AS 03-10-75// > 
ACCEPT 6, DATE!, DATE2 

6 FORMAT 2A4 ) 

C*****SPECIFY FLIGHT CONDITION FOR ENG INF 

7 TVF'E 11 


ACCEPT 12, F'O, TO, MO 
C :<:+ + +■ + INITIALIZE DACS 
1O0 CAM. OS CL O, I ERR) 

CALL QSCM, I ERR) 

105 CALL QWBDAS < V 1 , 0, 24, I ERR ) 
CALL. QSTDA 
DO 106 K=l, 16 


K"K=K— 1 

CALL QWCLLCKK, . FALSE , I ERR) 
106 CONTINUE 


11 FORMAT </3X, 33HTVPE DESIRED VALUES FOR P0, TO, M0 /) 

12 FORMAT <F7. 3, FS. 3, F5. 2) 

F'OS I = P0*. 6394 S 

T 0S I = T0+. 55555 
V3= P0/20. 

V5-M0/3. 


V6=T 8/1000, 

V 1 1 - V5+ S S Q R T < VS ) 

C*****PLACE analog in ic mode 

19 TVF'E 20 

20 FORMAT OX, 48HSLAVE CONSOLE 2 TO CONSOLE 1. MANUALLY GO TO If: /) 
TVF'E 21 

21 FORMAT C/3X, 44HPR0CEED TO DYNAMIC PART OF PROGRAM BY R-S-R /) 
PAUSE 

C*****READ ADC VALUES AND GENERATE MAP OUTPUTS 

22 CALL GRBADSCXG, 0, 5, I ERR) 

Y7=MAF'2C4, X3, X4> 



V5=MAP2L<7) 

IFCX2. GT. ©S> X2=. 0S 


A=MRP2vl, X4, -X2 ) 

V7=<Y7*<. 5S+. 5S*A) )/. 5b 
DV7=. 0S 

REV= < . 56715S*X0*<. 9S*X1+. 
I F < REV. LT . 30 5 84 S > DV7=~. 
I F < REV. LT. 14202S> DV7=- 
V7=V7-DV7 


17874S ) ) / <. XI* XL ) 
02626S*REV+. 00803S 
15077S+REV+. 02571S 


V3=V7 

CALL QUJDAS<V3, 3, IERR> 

CALL QWJDASCVS, 5, I ERR > 

CALL QUJDAS<Y7, 7 , I ERR ) 

CAI.L OREADS < X5, 5, 3, I ERR ) 

IF*' X5. GT, 00008S> X 5= 00000S 


V4=MAP2 < 2.. X6, 


V€=MAP2C 


{ .6 ) 


Y6«<V6*<. 5S+. 5S+V4 ) ) /. 5b 


V2=V6 

CALL QWJDRSCV2, 2, I ERR) 
CALL QUJDAS<Y6, 6, IERR) 


223 CALL QRBADSCXS, 8- 3- I ERR) 


RT4=SSQRT<X8) 

V7= < . 7 956 7 S* XI 0 ) XRT 4 
V8-MAP2<5, X9, V7) 
V9=MAP2L<8) 

V9= < -V9*RT4 > X. 78259b 
CALL QWJDAS<Y8, 8, I ERR) 


CALL QWJBASCV9, 9, I ERR) 

224 CALL QRBADSOili, 11, 3, I ERR) 
RT41-SS0RT < X12 ) 

VS=<. 8 2 6 9 5 S + X 1 1 ) X R T 4 i 
V10=MAP2<S, XI 3/ V8) 

V10=MRP2L<9) 

Vll= <■ — V10*RT 41 ) /. 46715b 
CALL QWJDASCVIG, 10.. IERR) 

CALL OUJDASvVil, 11, IERR) 

225 CALL OREADS <X22, 22, 2, IERR) 
V4=X22*X23 

CALL QWJDASCV4, 4, IERR) 

226 CALL ORBADS<X14, 14. 5, IERR) 

IF<V5. GT. . 36667S ) GO TO 2265 

AE=. 04162S+X14-. 08124S+. 73465S+X17 
GO TO 2266 


2265 

2266 


AE=. 0 25 8 5 S * X 14 — . 05837S+. 80875S+X17 
AR=AE/X17 

DPR-C 69683S*AR+AR-AR+. 32094b)/. 10567b 
PRSUB= . 52828S-DPR 


PE-PRSUB+X18 

FRD=. 7 5 7 2 2 S * X 16+ V 1 1 

IFC. 20000S+V3. L-T. PE) GO TO 227 

PRE=. 2 0 0 0 0 S * V 3 / X 1 8 

V12-<! <X:1 5*SSQRT < i. . 2578S- 2578S*PRE)+X14) )/. 
1 /. 34933S 


7652S-FRD) 


GO TO 228 


227 PRSUP=. 52828S+DPR 
PE=PRSUP*X18 

VI 2= < <X15*SSGRT< <. 47302S-PRSUP/. 90106S+. 
1 •+ X14 ) ) /. 69633S— FRD— < AE+ < . 20000S*V3-PE ) ) 


85065S+PRSUP + PRSUP 
j:- 4 y s lz> 
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228 V7SQ=V7+V? 

PRF1= 76979 S + V 7 SQ + . 01157S+V7+. 13806S 
PRF2=<. 9 G 6 4 5 S + V 7 S Q - , 41427S+V7+. 18519S)/ 8S 
V0=C 75 S+PRF ± + < 99999S+X2 >— X2+PRF2 ) /. 75S 
V0~ V0 — . 83333S+X3 > /V0 
V6SQ=V6+V6 

PRCi= < . 291 1. 7 S + V6S Q- . 46385S+. 7167S+Y6)/. 65 
PRC2=<. 8S275S+V6SG!—. 80752S+Y6+. 245495)/. 25S 
VI = 6S+PRC1+C. 99999 5 + X 5 > - X 5 + P R C 2 > / . 6S 
VI = < VI- 83333S+X7 ) /VI 
22. CALL QWJDASCV12, 12, I ERR) 

CALL QW JDAS < V0, 0, I ERR) 

CALL QW JDAS < VI, 1, I ERR > 

C + + + + + 0 U J P U T UNSCALED DATA AT TELETYPE IF DESIRED 
IF<. NOT SENSW < 1 ) ) GO TO 22 
CALL QRAMICILOC) 

CALL QSC < 2, IERR> 

CALL QSH < I ERR ) 

CALL QSC<0, I ERR> 

CAI.L QSC Cl, I ERR) 

CALL QRBA0SCX19, 19, 2, I ERR) 

CALL QSC C 2, I ERR) 

P13Q2-X3 

P13G2-P13G2+1. 5+3. 1064 
XNL.R2=X4 

XNLR2-"XNLR2+i. 2+10289. 

P3G22-X7 

P3Q22=P3Q22+1. 5+8. 4113 
XNHR22=X6 

X N H R 2 2 = X N H R 2 2 + 1 . 1+16777. 

P41Q4=X9 

P41Q4=P41Q4+3. 0+, 27448 
CNHPT = V7 

CNHPT =CNHPT +1. 25+238. 47 
CHHTSI -CNHPT +1. 3418 
P5Q41=X13 

P 5 Q4 1 = P5Q4 1 +2 . 5+. 44355 
CNLPT = V8 

CNLPT -CNLPT +1. 5+220. 80 
CNLTSI = CNLPT +1. 3416 
GVIP0S-X2 
GV I P03=GV I PCS +25. 

H7SP0S-X5 

H VSP05= H VSPOS+4 4 . +4. 0 

P2-X0 

P2-P2+40. 

P2SI =F’2+. 68948 
T2~X1 

T2-T2+10G0. 

T2SI=T2+. 55555 
NF4-X22 

NF4-NF4+4. 5833 
WF4SI = WF4+. 45359 
WF4=WF4+3600. 

AN" XI 7 

AN = AN +1000. 

ANSI ~ AN+. 00664516 
AN- AN/144. 



P4--P4 + 600 
F*4S I =■•• P4 + 6894::.: 

HNL^H.U 

: :hi r- hni +.15000, 

XNK»X10 

: :nh- xnh+15000 
T > i -xs 

T4“- T4 + 4060. 

T4SI- T4+. 55555 
PLR - HI 9 
pi n- PI.R+-J50. 

T7~H14 

T7 = T 7+50Q0 
T 7 S I -- T 7 + 5 5 5 5 5 
MF7 = K26 
MF7 -• WF7+20. 

WF7SI= WF7*. 45253 
WF 7 ~ WF7 : + 2 SO0 . 

WA?=H16 
M Fi r - W Fl 2 + 4 5 G 
UA2SI--WA2+. 45253 
T 41™ HI 2 
T 4 1= T 4 1+2000. 

T41SI “*T4i+. 55555 
P21Q2- V5 

P 2 1 Q 2 == P 2 1 0 2 •+ 1 . 5 + 2. 0986 
WAR22™ Y6 

WAR?2=WRR22+1 1+55. 697 
WR22S I - WAR 22+ 45259 
WAR2 == ; V7 

WAR 2 ~ N fl R 2+ 1 . 2+220. 28 
WAR2S I = WAR 2*. 45253 
L-JPHPT” VS 

WPHF’T -WF’HPT +1. 50+ . 072872 
WF'HTSI = WPHPT : +. 26548 
HP HPT * V3 

HF'HPT =HPHPT +1. 50+. 28660 
HPHTSI* HPHF’T+2118. 7 
WPLPT= Vi© 

WPLPT-WF‘LF'T+2 5*. 29226 
WPLTS I =WF*LPT +. 26548 
HPLPT“ ; V10 

HPL PT=HPLPT +2. 5 +. 29852 
HALTS I = HPLPT+2118. 7 
FN-V12 
FN~ FN-+ 2 000 0 . 

FNS 1 = r Fh4+4 . 4482E-3 


REV I =REY 


F'.EV I =REV I +2. 

TV PE 24/ DATE!.. DATE2 

24 FORMAT 21H/ 22HF106 SIMULATION 
1 2R4/V/0 

IFC NOT. SENSW<2>> GO TO 441 
TYPE 25.- P0SI . P0 

25 FORMAT < 5X, 9HP0 = , F7, 2/ 9X, 
1 > 


STEADY-STATE DATA.. 5K/ 5HDATE : 


7 HN/SQ CM/ 7H/ 2Hv F7. 2.* 9H.* SHP£ 


TYPE 26/T0SI..T0 
26 FORMAT 5H/ 3HT0 


/ F7 


7 HR 


.. 7H.« 2H( ■ F7. 2/ 9X, 8 HR 


r- ir 


28 



« ox :?hi so ■ • r .7 * 

i *iT- f .0 , P2'"I iv 

•V FORMAT • OX, ?i |py • f 7. 3X, 7HN/50 CM, ?X, 2H< , F7. 3, 9X, 8HPSIA > 

t s f t. ivi , rx 

; - F OP MAT ■: OX 9UT ■’ - , FT 2, 3X, 7HK . 7X, 2Hv , F7. 2, 9X, 8 HR > 


TYPE 2 : <0, F'FVi 


•■0 F OR HR T «: V’, 3HPEYI 

, 

FT 

4 > 


TS’F'E CO.Fl.R 
70 FORMAT 9HF- LF: 


F7 

2.* 9X, 

3HDEG > 

TVPF 71 . WF4SI . Mr 4 
A r c^ifa , 9HHF4 

1* 

TYPE S2. MF7SI . MF7 
72 FOR HR 7 3 OX, Of IMF 7 

- .■ 

F7. 

4, 9X, 

7 HI:. GX SEC , 7X, 2H< 


F7. 

,- A , v , 

7HKG/SEC , 7X • 2FK 

i.» 

TV HR GVIPOS 

22 FOPl’iFlT ' OX ■ 9HGV I PCS 


FT, 

3: , 9X • 

SHDEG;- 

! V FF 24 . HVSPOS 
74 FORMAT', OX, 9HHVSP0S 


F7 

2 t Q-, _ 

3hdeg:> 

TVPE 373, ANSI ■ RN 
3 0 FORMAT f OX, 3 HAN 

-• 

FT. 

5, 9X, 

7HSQ H , 7X, 2HC 

17 

TVPE 26, XNH 
36 FORliRT '’OX, 9HXNH 


F7. 

0, 9X, 

3HRFT1 > 

TVPE 37-XNL 
37 FORMAT ‘OX, 9HXNL 

r- 

F7. 

0, 9X, 

3HRF1T > 

TVPE 28, WH2SI, MR? 
3 8 FORMAT >:! OX, 9HWR2 

=r. 

F7. 

2, 9X, 

7HKGXSEC , 7X, 2H 

1 ‘ > 

TVPE 40 i P4SI, P4 
■ 1 0 FORMA T*.’ OX, 9 HP 4 


F7 

Q* * 

7HN/SQ CM, 7X, 2H< 

I 

TVPE 41, FNSI-FN 
41 FORMAT < OX, 9HFN 


F7. 

2 j 9X, 

7HKN , 7X, 2H< 

1 > 

TVPE 42, T4SI .T4 
4 2 FORMAT < 5 X , 9 H T 4 


FT. 

1, 9X, 

7HK , 7X, 2HC 

1 > 

TVPE 39, T41SI, T41. 
39 FORMAT FOX, 9HT41 

= , 

F7. 

1, 9X, 

7HK , 7X, 2HC 

TVPE 4 3:, T7SI,T7 
4 7 FORMAT COX, 9HT7 


F7. 

1, 9X, 

7HK , 7X, 2H< 

1 ■ 

TVPE 40, Pi 202 
I PC NOT SENSPK3’; 1 > 

GO 

TO 

60 


40 FORMAT C 5! 7 9 HP 13Q2 

= , 

F7. 

4 > 


41 TVPF 46.XNRR2 
4 6 FORMAT C OX, 9HXNLR2 

XU 

F7. 

0, 9X, 

3HRPM > 

TVPE 47, WAR2SI, WAR 
47 FORMAT C OX, 9HWAR2 

2 

F7 

2, 9X, 

7HKG/SEC , 7X, 2HC 


1 > 

TVPE 48, F21Q2 


, F7. 8, 9X, 8HLBM/HR > 
, F7 0, 9X, 8HLBM/HR > 


F7. 4, 9X, SHSQ FT 


, F7. 2, 9X, 8HLBM/SEC ) 
, F7. 2, 9X, SHF'S I fl > 


, F7. 0, 9X, 8HLBF > 

, F7. 1, 9X, 8 HR > 

, F7. 1, 9X, 8HR > 

, F7. 1, 9X, 8 HR > 


, F7. 2, 9X, 8HLBM/SEC ’> 


29 



48 FORMAT < 5X, 9HP21.Q2 = 

TVPE 49, P3Q22 

49 FORMAT CSX, 9HP3Q22 = 

TVPE 56, XNHR22 

56 FORMAT CSX, 9HXNHR22 = 
TVPE 51, WR22SI, WAR22 

51 FORMAT CSX, 9HWAR22 * 
1> 

TVPE 52, P41Q4 

52 FORMAT < 5X, 9HP41Q4 = 

TVPE 53, CNHTSI, CNHPT 

53 FORMAT <5X, 9HCNHPT = 
1> 

TVPE 54, WPHTSI, WPHPT 

54 FORMAT C 5X, 9HWFHPT = 
1 > 

TVPE 55, HPHTSI, HP HPT 

55 FORMAT CSX, 9HHPHPT = 
1> 

TVPE 56, P5Q41 

56 FORMAT C 5X, 9HP5Q41 - 

TVPE 57, CNLTSI, CHLPT 

57 FORMAT < 5X, 9HCNLPT - 
1> 

TVPE 58, WF'LTSI , UP LPT 

58 FORMAT <5X, 9HWPLPT 
1> 

TVPE 59, HPLTSI, HPLPT 

59 FORMAT < 5X, 9HHPLPT = 
1> 

TVPE 591. V0 

591 FORMAT <5X,' 9HSMF 
TVPE 592, VI 

592 FORMAT C 5X, 9HSMC 

60 IFC I LOG. EQ. 6 > GO TO 
CALL Q50PC I ERR) 

GO TO 62 

61 CALL QSICCIERR> 

62 CALL QSCv 0, I ERR > 

CALL QSC < 1 ■ I ERR > 

GO TO 22 

END 


, F7. 4 
, F7. 3 > 

, F7. 0, 9X, 3HRPH > 

. F7. 2, 9X, 7HKG/SEC , 7X, 2HC 

, F7. 5> 

, F7. 2, 9X, 7H , 7X, 2HC 

, F7. 5, 9X, 7H , 7X, 2HC 

, F 7. 1, 9X, 7H , 7X, 2H< 

, F7. 5> 

, F7. 2, 9X, 7H , 7X, 2HC 

.. F7. 5, 9X, 7H , 7X, 2HC 

, F7 1, 9X, 7H , 7X, 2HC 


61 


F7. 2, 


SHLBM/5EO 


, F7 2, 9X, 8H 
, F7. 5, 3X, 8H 
, F7. 5, 9X, 8H 


, F7. 2, 9X, SH 
, F7. 5, 9X. 8H 
, F7. 5, 9X, 8H 
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FORTRAN Symbols 


A 

AE 

AN 

ANSI 

AR 

CNHPT 

CNHTSI 

CNLPT 

CNLTSI 

DPR 

DY7 

FN 

FNSI 

FRD 

GVIPOS 

HP HPT 

HPHTSI 

HPLPT 

HPLTSI 

HVSPOS 

IERR 

DC 

JY 

K 

KK 

KX 

KY 

MAP2 


shift in fan map corrected airflow due to change in inlet guide vane position 
(scaled) 

exhaust nozzle exit area (scaled) 

2 

exhaust nozzle throat area, ft 

2 

exhaust nozzle throat area, m 

exhaust nozzle expansion ratio (scaled) 

high-pressure-turbine corrected speed, 

high-pressure- turbine corrected speed, 

low-pressure-turbine corrected speed, rpm/ u R i/<s 

1/2 

low-pressure-turbine corrected speed, rpm/K ' 
shift in critical pressure ratio due to expansion ratio 

shift in fan map corrected airflow due to change in Reynolds number (scaled) 

net thrust (uninstalled), lbf 

net thrust (uninstalled), kN 

ram drag (scaled) 

inlet guide vane position, deg 

high-pressure-turbine enthalpy drop parameter, Btu/lbm - °R * / 2 - r pm 
high-pressure-turbine enthalpy drop parameter, J/kg-K ' -rpm 
low-pressure-turbine enthalpy drop parameter, Btu/lbm- °R*/ 2 -rpm 
low-pressure-turbine enthalpy drop parameter, J/kg-K ' -rpm 
stator vane position, deg 
error flag for linkage routines 

array containing number of points per curve for each map pair 
array containing number of curves for each map pair 
control line initialization index 
K-l 

array containing x out- of- range counts for each map pair 
array containing y out-of-range counts for each map pair 
bivariate function (first function) 


rpm/°R 1//2 

rpm/K */ 2 


31 



MAP2L 

MO 

NX 

NY 

PE 

PLA 

PRC1 

PRC2 

PRE 

PR FI 

PRF2 

PR SUB 

PRSUP 

P I 

PJSI 

PJQI 

QRAMI 

QRBADS 

QPC 

QSH 

QSIC 

QSHYIN 

QSOP 

QSTDA 

QWBDAS 

QWCLL 

QWJDAS 

REY 

REYI 

RT4 

32 


bivariate function (second function) 

Mach number 

array containing nu tiber of points per curve for each map pair 
array containing number of curves for each map pair 
exhaust plane pressure (scaled) 
power lever angle, deg 

compressor critical pressure ratio, axial vanes (scaled) 
compressor critical pressure ratio, cambered vanes (scaled) 
nozzle pressure ratio 

fan critical pressure ratio, axial vanes (scaled) 
fan critical pressure ratio, cambered vanes (scaled) 
critical nozzle pressure ratio 
design pressure ratio for supersonic nozzle flow 
pressure at station I, psia 

2 

pressure at station I, N/cm 

ratio of pressure at station J to pressure at station I 

linkage routine for sensing analog mode 

linkage routine for reading ADC's 

linkage routine for selecting analog console 

linkage routine for placing analog console in HOLD mode 

linkage routine for placing analog console in IC mode 

linkage routine for addressing analog consoles 

linkage routine for placing analog console in OPERATE mode 

linkage routine for transferring DAC data 

linkage routine for loading DAC's 

linkage routine for setting control lines 

linkage routine for "JAMMING" DAC's 

Reynold's number index (scaled) 

Reynolds number index 
square root of T4 (scaled) 



RT41 

SENSW 

SSQRT 

T I 

TKI 

TJQI 

V3 

V4 

V5 

V6 

V7 

V8 

V9 

V10 

VI 1 

WA2 

WA2SI 

WAR2 

WAR22 

WAR2SI 

WF4 

WF4SI 

WF7 

WF7SI 

WPHPT 

WPHTSI 

WPLPT 

WPLTSI 

WR22SI 


square root of T41 (scaled) 

array containing logical indication of sense switch positions 
scaled-fraction square root routine 
temperature at station I, °R 
temperature of station I, K 

ratio of temperature at station J to temperature at station I 
ambient pressure (scaled) 

shift in compressor map corrected airflow due to change in stator vane posi- 
tion (scaled) 

Mach number (scaled) 

ambient temperature (scaled) 

high-pressure-turbine corrected speed (scaled) 

low- pres sure- turbine corrected speed (scaled) 

high-pressure-turbine enthalpy drop parameter (scaled) 

low- pres sure- turbine enthalpy drop parameter (scaled) 

product of Mach number and square root of ambient temperature (scaled) 

fan airflow, lbm/sec 

fan airflow, kg/sec 

fan corrected airflow, lbm/sec 

compressor corrected airflow, lbm/sec 

fan corrected airflow, kg/sec 

main- combustor fuel flow, lbm/hr 

main- combustor fuel flow, kg/hr 

augmentor fuel flow, lbm/hr 

augmentor fuel flow, kg/hr 

high-pressure-turbine corrected flow, lbm-°R-in. Vlbf-rpm-sec 
high-pressure- turbine corrected flow, kg-K-cm /N-rpm-sec 

n 

low-pressure-turbine corrected flow, lbm-R-in. /lbf-rpm-sec 
low-pressure-turbine corrected flow, kg-K-cm /N-rpm-sec 
compressor corrected airflow, kg/sec 
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XNH 

high-speed- rotor speed, rpm 


XNHR22 

compressor corrected speed, rpm 


XNL 

low- speed-rotor speed, rpm 


XNLR2 

fan corrected speed, rpm 


XVALS 

array containing scaled map input x 

data 

XI 

variable read on ADC channel I 


YVALS 

array containing scaled map input y 

data 

YJ 

variable output of DAC channel I 


Y6SQ 

output of DAC channel 6 squared 


Y7SQ 

output of DAC channel 7 squared 


ZVALS 

array containing scaled map output z 

data 
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Fan - console 1 





Compressor - console ! 1 
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High pressure turbine and inter-turbine volume - console 2 





Low pressure turtiine and mixing volume - console 2 



45 000 
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Augmentor - console 2 


(AH>lt^ w 4 1 + 0. 01436 KblwlT^ 2 ) 



Torque balance - console 2 
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Miscellaneous - console 2 



REFERENCES 


1. Stone, C. R.; Miller, N. E. ; Ward, M. D. ; and Schmidt, R. D. : Turbine Engine 

Control Synthesis. Vol. I: Optimal Controller Synthesis and Demonstration Final 
Technical Report. AFAPL-TR-75- 14- VOL- 1, Honeywell, Inc. (AD-A0142 2 9), 
1975. 

2. Michael, Gerald J. ; and Farrar, Florence A. : An Analytical Method for Synthesis 

of Nonlinear Multivariable Feedback Control. UARL-M941338-2, United Aircraft 
Corp. (AD-762797), 1973. 

3. Michael, Gerald J. ; and Farrar, Florence A. : Development of Optimal Control 

Modes for Advanced Technology Propulsion Systems. UARL-M911620- 1, United 
Aircraft Corp. (AD-767425), 1974. 

4. Bowles, Robert J. : Sub-Optimal Control of a Gas Turbine Engine. GE/EE/73A-1, 

Air Force Inst. Tech. (AD- 777852), 1973. 

5. Bentz, C. E. ; and Zeller, J. R.: Integrated Propulsion Control System Program. 

SAE Paper 730359, April 1973. 

6. Szuch, John R.; and Bruton, William M. : Real-Time Simulation of the TF30-P-3 

Turbofan Engine Using A Hybrid Computer. NASA TM X-3106, 1974. 

7. Szuch, JohnR.; and Seldner, Kurt: Real-Time Simulation of FI 00- PW- 100 Turbo- 

fan Engine Using the Hybrid Computer. NASA TM X-3261, 1975. 

8. Integrated Propulsion Control System (IPCS). Volume H - Technical Description. 

Boeing Aerospace Co. (AFAPL-TR-76-61), 1976. 

9. Cwynar, David S . ; and Batterton, Peter G. : Digital Implementation of the 

TF30-P-3 Turbofan Engine Control. NASA TM X-3105, 1974. 

10. Deiters, R. M. ; and Nomura, T. : Circle Test Evaluation of Method of Compensat- 

ing Hybrid Computing Error by Predicted Integral. Simulation, vol. 8, no. 1, 
Jan. 1967, pp. 33-40. 

11. Bekey, George A. ; and Karplus, Walter J. : Hybrid Computation. John Wiley & 

Sons, Inc., 1968, pp. 117-124. 

12. Szuch, John R. : Application of Real-Time Engine Simulations to the Development 

of Propulsion System Controls. NASA TM X-71764, 1975. 


46 



TABLE I. - DESIGN PARAMETERS 


Compressor discharge volume, V 3 , m 3 

0. 0468 

Main- combustor volume, V 4 , m 3 

0. 0468 

Interturbine volume, V 4 m 3 

0. 6561 

Mixing volume, Vg, m 3 

0.8470 

Augmentor volume, V 7 , m 3 

0.7128 

Duct volume, V^ 3 , m 3 

1.427 

Augmentor inductance, (*/Ag c ) AB , N-sec 2 /kg-cm 2 

0.0007598 

Duct inductance, (f/Ag c ) D , N-sec 2 /kg-cm 2 

0. 0007598 

ffigh-speed rotor inertia, I H , N-cm-sec 2 

565.35 

Low-speed rotor inertia, I L , N-cm-sec 2 

610.00 

Main- combustor pressure loss coefficient, Kg, N 2 sec 2 /cm 4 -K-kg 2 

0.00114 

Low-pressure- turbine discharge pressure loss coefficient, K pR5 

1. 024 

Augmentor pressure loss coefficient, K AB , N 2 -sec 2 /cm 4 -K-kg2 

3. 5659xlCf 6 

Nozzle flow coefficient, K^, kg-K 1//2 /N-sec 

0. 1509 

Fraction of high-pressure- turbine cooling bleed that performs work, K DX „, T ™ 

0.6292 

Fraction of low -pres sure -turbine cooling bleed that performs work, K B ^ WLr p 

0. 1114 

Fan inlet specific heat, c p 2 , J/kg-K 

1009 

Compressor inlet specific heat, c p 2 2 , JAg-K 

1001 

Compressor discharge specific heat, c p 3 , JAg-K 

1039 

Main- combustor specific heat, c p 4 , J/kg-K 

1145 

Interturbine specific heat, c p 4 p JAg-K 

1116 

Mixing-volume specific heat, c p g, J/kg-K 

1062 

Augmentor specific heat, c p 7 , J/kg-K 

1062 

Duct inlet specific heat, c p 13 „ JAg-K 

1009 

Duct discharge specific heat, c p 16 , J/kg-K 

1030 

Main- combustor specific heat ratio, y 4 

a l . 292 

Interturbine specific heat ratio, y 4 1 

1.306 

Mixing volume specific heat ratio, yg 

1.344 

Augmentor specific heat ratio, y 7 

b l. 359 

Compressor discharge temperature time constant, r^, sec 

0.05 

Duct temperature time constant, r^, sec 

0.05 

Heating value, HVF, J/kg 

4.407xl0 7 


Effectively decreased by a factor of 20 to match baseline digital data. 
Effectively decreased by a factor of 10 to increase simulation stability. 
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Inlet 




Figure 2. - Computational flow diagram of real-time FI 00- PW- 100(3) engine simulation. 





















Fraction of design compressor 
corrected speed 
1.0672 




Figure 6. - Effect of variable stator vane position on F- 100- PW- 100(3) compressor performance 
map. 






Fraction of design low-pressure-turbine enthalpy Fraction of design low-pressure-turbine flow parameter 

parameter, fo 



Fraction of design low-pressure-turbine pressure ratio 
(b) Enthalpy parameter. 

Figure 8. - F100-PW- 100(3) iow-pressure-turbine performance maps. 


Augmentor fuel-air ratio 

Figure 9. - F100-PW- 100(3) augmentor effi- 
ciency function. 



Fraction of design duct airflow (corrected to station 16) 
Figure 10. - FIOO-PW- 100(3) duct pressure loss function 



















Fraction of design main combustor temperature 













.6 .7 .8 .9 1.0 1.1 .8 .9 1.0 


Fraction of design compressor corrected airflow 

(c) Altitude, 9. 144 kilometers; Mach number, 0.9. (d) Altitude, 13. 72 kilometers; Mach number. 0.9. 

Figure 19. - Comparison of open-loop hybrid and baseline digital steady-state data for compressor operating line at standard-day conditions. 
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